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摘 要 : 全 球 变 暖 导致 干旱 等 极端 事件 频 发 ,严重 威胁 生态 安全 和 社会 经 济 可 持续 发 展 , 尤 其 是 于 
旱 区 这 一 气候 响应 敏感 区 。 本 研究 基于 流域 站 点 观测 气象 ,水 文 等 数据 和 第 六 次 国际 耦合 模式 比 
较 计 划 (CMIP6) 全 球 气候 模式 数据 ,应 用 分 布 式 水 文 模型 HEC-HMS 模型 ,模拟 并 预 估 了 塔里木 河 
重要 源流 一 一 叶 尔 羌 河流 域 历 史 (1986 一 2014 年 ) 与 未 来 (2015 一 2100 年 ) 径 流 变 化 趋势 以 及 水 文 


干旱 风险 。 结 果 表 明 :(1) HEC-HMS 模型 在 干旱 区 流域 适用 性 较 好 ,未 来 3 个 共享 社会 经 济 路 径 
(SSPs) 下 叶 尔 羌 河流 域 年 径流 和 标准 化 径流 指数 (SRI) 都 将 呈 显 著 增 加 趋势 (P<0.1),SRI LAR 
3€ 2) 0.13-0.27 - (103), (2) 构建 并 对 比 了 流域 历史 和 未 来 时 期 4 个 干旱 特征 变量 的 边缘 分 布 , 未 来 
干旱 持续 时 间 干旱 强度 和 烈度 峰值 都 明显 大 于 历史 时 期 ,干旱 程度 可 能 将 不 断 加 剧 ,(3) 相 较 于 


历史 时 期 , 叶 尔 羌 河 流域 未 来 水 文 干旱 发 生 的 联合 概率 将 持续 减 小 , 且 未 来 干旱 事件 的 重 现 期 都 
将 有 不 同 程度 的 延长 。 研 究 结 果 可 为 叶 尔 羌 河 流域 水 资源 管理 与 规避 水 文 干 旱 风 险 措施 的 制定 


提供 科学 参考 。 
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全 球 气候 变化 加 速 了 水 循环 ,导致 包括 干旱 在 
内 的 极端 水 文 事件 频 发 ,尤其 是 北半球 中 纬度 地 
区 “与 此 同时 ,人 口 增长 和 农业 、 能 源 和 工业 部 
门 的 扩张 等 其 他 因素 ,预计 将 进一步 加 剧 水 资源 短 
缺 ** ,特别 是 在 水 资源 匮乏 的 干旱 地 区 。 水 文 干旱 
是 指 由 于 降水 的 长 期 短缺 而 造成 某 段 时 间 内 地 表 
水 或 地 下 水 收文 不 平衡 ,出 现 水 分 短缺 ,使 江河 流 
量 湖泊 水 位 、 水 库 蓄 水 等 减少 的 现象 ,其 发 生 会 对 
区 域 的 农业 、 生 态 和 社会 经 济 发 展 产 生 重要 的 影 
响 。 相 较 于 农业 、 气 象 等 干旱 ,水文 干旱 特征 更 加 
复杂 ,所 产生 的 后 果 也 更 为 严重 ""。 河 流 径 流 是 干 
早 区 最 主要 的 水 资源 ,维持 着 下 游 绿洲 农业 和 荒漠 
生态 系统 的 发 展 ,水 文 干旱 能 更 全 面 地 反映 以 灌溉 
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农业 为 主 的 干旱 区 的 干旱 情况 ,因此 更 需要 在 对 未 
来 径流 过 程 进行 可 靠 模拟 与 预测 的 基础 上 ,进行 流 
域 未 来 水 文 十 旱情 况 分 析 和 预 佑 ,以 期 为 地 区 水 
资源 有 效 的 规划 和 管理 提供 科学 有 效 的 依据 。 

内 外 学 者 对 水 文 干旱 做 了 大 量 研 究 , 包 括 水 
文 干旱 指标 的 构建 与 特征 变量 的 识别 .水 文 干旱 的 
形成 机 制 .历史 与 未 来 水 文 干旱 风险 评估 等 ,研究 
尺度 包括 流域 尺度 到 全 球 尺 度 ,技术 方法 日 趋 成 
熟 。 未 来 径流 的 模拟 和 预测 是 未 来 水 文 干 旱 风 险 
分 析 的 基础 ,而 水 文 模型 是 研究 流域 径流 过 程 的 重 
HTH, HEC-HMS 是 一 种 具有 一 定 物 理 意义 
的 分 布 式 水 文 模 型 , 旨 在 模拟 流域 尺度 的 降水 - 融 
雪 径 流 过 程 ,已 被 广泛 使 用 。 例 如 ,Azmat RK 
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模型 应 用 于 喜马拉雅 山脉 西部 由 冰雪 融 水 补给 的 
高 海拔 Jhelum 河 流域 ,并 预测 了 未 来 气候 变化 情景 
下 的 径流 量变 化 。 未 来 径流 的 预测 主要 基于 率 定 
好 的 水 文 模型 ,通过 全 球 气 候 模 式 获 得 的 未 来 气候 
变化 情景 数据 作为 输入 "”"""。2012 年 秋 , 由 世界 气 
候 研 究 项 目 “ 气 候 模 拟 工 作 组 ”组 织 的 全 球 耦 合 模 
式 比较 计划 第 六 阶段 (CMIP6) 正 式 启 动 , 以 共享 社 
会 经 济 路 径 (SSPs ) 为 核心 的 新 社会 经 济 气候 变化 情 
景 反映 了 辐射 强迫 与 社会 经 济 发 展 的 关系 ,将 成 为 
IPCC 最 新 报告 的 重要 理论 基础 ”。 

叶 尔 羌 河 是 塔里木 河 最 长 的 支流 ,流域 内 全 年 
降水 稀少 ,下 游荡 漠 生 态 系统 与 绿洲 农业 用 水 全 部 
依赖 于 上 游 来 水 ,上 游 径流 补给 来 源 主要 为 冰雪 融 
水 和 山区 降水 ,径流 年 际 和 年 内 变化 大 。 受 这 一 特 
殊 的 补给 方式 影响 ,流域 水 文 过 程 对 气候 变化 极为 
敏感 ,上 游 的 水 文 格局 对 下 游 及 塔里木 河 干流 的 来 
水 和 用 水 情况 有 很 大 的 影响 ,而 水 文 干旱 对 该 区 域 
的 影响 非常 严重 ,尤其 以 3 一 5 月 的 春 旱 最 为 严重 。 
随 着 下 游人 口 .耕地 等 的 扩张 ,水 资源 需求 更 大 ,与 
此 同时 ,气候 变化 引起 的 极端 气候 水 文 事件 增多 ， 
势必 会 影响 塔里木 河流 域 水 资源 供需 平衡 ,威胁 整 
个 塔里木 河流 域 水 安全 。 目 前 ,已 有 学 者 对 叶 尔 羌 
河流 域 历史 时 期 水 文 过 程 及 水 文 干旱 特征 进行 了 
研究 ,但 缺乏 对 未 来 时 期 水 文 干旱 概率 及 风险 的 预 
估 。 本 研究 以 塔里木 河 重 要 源流 一 一 叶 尔 芜 河 为 
研究 区 ,应 用 分 布 式 水 文 模型 HEC-HMS 模 型 ,结合 
6 个 CMIP6 全 球 气候 模式 ,选择 其 中 的 3 个 SSPs 情 
景 ,估算 了 流域 控制 期 (1986 一 2014 年 ) 和 未 来 
(2015 一 2100 年 ) 径 流 数据 ,进一步 计算 标准 化 径流 
指数 (SRI) 作 为 未 来 水 文王 旱 评 价 指标 ,借助 游程 
理论 和 多 维 Copula 函数 , 预 估 叶 尔 羌 河流 域 在 未 来 
不 同 SSPs 情景 下 各 个 水 文 干 旱 特征 变量 的 变化 趋 
势 ,以 及 不 同情 景 下 联合 干旱 概率 以 及 重 现 期 ,为 
政府 和 相关 部 门 应 对 气候 变化 ,制定 相应 的 水 资源 
管理 措施 ,更 好 地 服务 于 “一 带 一 路 ”建设 提供 科学 
依据 。 


1 数据 与 方法 


1.1 研究 区 概况 
叶 尔 羌 河 流域 (76°40'~79°25'E,37°20'~40°10'N) 
位 于 塔里木 河流 域 西南 ,发 源 于 喀 喇 昆仑 山 ,流域 
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图 1 研究 区 示意 图 
Fig. 1 Schematic diagram of the study area 


出 山口 以 上 集 水 面积 46704 km? ,平均 海拔 4400 m, 
是 塔里木 河 最 长 的 源流 22 (图 1)。 叶 尔 羌 河 流域 
位 于 典型 的 内 陆 干旱 气候 区 ,在 河流 的 上 游 山区 发 
育 着 现代 冰川 和 永久 积 雪 。 出 山口 以 上 流域 内 冰 
川 面积 5852.57 km’, 占 流域 面积 的 11.95%, 积 雪 面 
积 约 占 流域 面积 的 15%”" 。 流 域 主要 由 冰雪 融 水 
补给 ,径流 量 年 际 变化 较 大 , 且 受 气候 变化 影响 显 
著 下 。 叶 尔 羌 河 下 游 灌 区 是 新 疆 最 大 的 农业 灌区 ， 
是 新 疆 力 至 全 国 最 重要 的 优质 粮 、 棉 、 果 品 生产 基 
地 四 ,流域 上 游 来 水 在 径流 供给 维持 下 游 生态 安 
全 和 社会 经 济 发 展 等 方面 至 关 重 要 。 
1.2 数据 来 源 与 处 理 

本 研究 选取 的 历史 时 期 气象 数据 为 站 点 观测 
数据 ,数据 来 自 中 国 国 家 气象 中 心 的 中 国 地 面 国际 
交换 站 气象 资料 日 值 数 据 集 ,包括 温度 .降水 .相对 
湿度 等 ,用 于 水 文 模型 的 构建 和 率 定 ;历史 时 期 水 
文 数据 来 自 叶 尔 羌 河 流域 出 山口 水 文 站 卡 群 站 观 
测 的 日 径流 数据 。 

用 于 预测 的 气象 数据 来 自 CMIP6 的 6 个 全 球 气 
候 模式 (BCC-CSM2,CNRM-CM6-1,CNRM-ESM2-1， 
GFDL-ESM4,MPLESM1-2-HR,MRI-ESM2-0) 输 出 的 
日 尺度 温度 、 降 水 等 数据 ,用 于 预 估 研 究 区 未 来 气候 
和 径流 变化 。 人 气候 模式 数据 时 间 尺 度 为 1986 一 2100 
年 ,其 中 1986 一 2014 年 为 控制 试验 期 ,2015 一 2100 
年 为 预 佑 期 ,选取 了 SSP126、SSP245 和 SSP370 3 个 
SSPs。 数 据 空 间 分 辨 率 都 为 230 km, 为 了 减 小 CMIP6 
全 球 气候 模式 数据 的 不 确定 性 ,本 研究 选择 Gauss- 
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ian 分 布 映 射 方法 (DM) 和 局 部 强度 缩放 方法 (LOCI) 
分 别 对 温度 和 降水 进行 降 尺 度 处 理 和 偏差 校正 ,这 2 
种 校正 方法 已 被 证 明 能 很 好 地 应 用 于 该 流域 的 气候 
模式 数据 校正 后 。 需 要 指出 的 是 ,由 于 CMIP6 气 候 模 
式 数据 中 控制 实验 期 为 1986 一 2014 年 , 故 本 研究 选 
取 的 气象 和 水 文 观测 数据 时 间 序 列 为 1986 一 2014 
年 ,用 于 对 比 验证 气候 模式 数据 的 适用 性 。 
13 研究 方法 
L3. 水 文 模型 HEC-HMS 模 型 是 美国 陆军 工程 兵 
团 水 文 工 程 中 心 (HEC) 开 发 的 流域 水 文 模拟 系统 
(HMS) ,是 一 个 具有 物理 概念 的 分 布 式 水 文 模型 。 
模型 输入 数据 包括 温度 .降水 . 蒸 散 发 .地 下 水 和 士 
壤 含 水 量 等 ,模型 中 的 融雪 模块 采用 温度 指数 法 来 
估算 积 雪 的 积累 和 消融 。 通 过 对 模型 中 7 个 模块 共 
26 个 参数 进行 率 定 ,得 到 最 优 模拟 结果 ,详细 模拟 
过 程 及 参数 见 Xiang 等 的 的 研究 。 本 研究 选取 无 量 
纲 的 拟 合 优 度 系 数 纳什 系数 (NSE) 、 相 对 偏差 
(PBIAS) 和 决定 性 系数 (RR)3 个 指标 对 模拟 结果 进 
行 评 价 。 在 水 文 模拟 中 ,如 果 NSE>0.75 IPBIASI< 
10% ,可 认为 模拟 效果 非常 好 ”。 
1.3.2 水 文 干旱 特征 分 析 本 研究 选取 3 个 月 尺度 
的 SRI 作 为 表征 水 文 干旱 程度 的 指标 。 参 考 前 人 研 
究 ,根据 不 同 的 SRI 阅 值 ,可 以 将 水 文 干旱 划分 为 
4 个 等 级 ( 表 1) ,本 研究 将 干旱 阔 值 设置 为 -0.5。 

采用 游程 理论 识别 干旱 事件 的 特征 变量 ,包括 
水 文 干旱 持续 时 间 (Dd) .水文 干旱 强度 (Ds) ,水文 
干旱 烈度 (Di) 和 水 文 干旱 烈度 峰值 (Dp)™。 
1.3.8 水 文 干旱 风险 估算 方法 Copula 函数 最 初 由 
Sklar 提 出 并 使 用 飞 。 利 用 边缘 分 布 和 相关 函数 ， 
Copula 函数 模型 可 以 构建 多 维 联合 分 布 ”2。 采 用 
二 维 三维 对 称 型 Archimedean Copula 函数 分 析 干 
旱 特 征 变量 之 间 的 关系 ,构造 干旱 历时 .干旱 强度 、 
干旱 烈度 和 干旱 烈度 峰值 的 联合 分 布 ,将 其 作为 评 
佑 干旱 风险 的 指标 ( 表 2)。 采 用 均 方 根 误差 (RMSE) 
和 赤 池 信息 准则 (AIC) 作 为 评价 各 种 Copula 函数 拟 
合 结果 优 劣 的 指标 。 

在 水 文 分 析 中 ,常用 重 现 期 来 度量 风险 的 大 
小 。 重 现 期 是 指 某 随 机 变量 的 取 值 在 长 时 期 内 平 
SYA TA) A RU, BE X, Xo AXP HI eS 
干旱 事件 中 具有 相关 关系 的 特征 变量 序列 ,边缘 分 


表 1 水 文 干旱 等 级 划分 


Tab.1 Classification of hydrological drought events 


标准 化 径流 指数 (SRI) 干旱 等 级 
(1 -0.5] 轻 度 干旱 
(-15,-1] PETE 
(72, -1.5] 重度 干旱 

<2 极端 干旱 


表 2 三 维 阿 基 米 德 Copula 函数 
Tab. 2 Symmetric Archimedean Copula functions 


Copula PR 公式 参数 范围 
Gumbel Ex(-(Chnu)' +(-Inu,)’ + (7j) [1, 99) 
Frank -Jp 人 + (e^ - De™ 一 (e^ - 1)/(e” 一 1j) iE 0), (0, 1] 
Clayton (u^ +u, +u,- 2)" (0, œ) 


1E :073 Copula 函数 的 参数 ;wu ww 为 干旱 变量 的 边缘 分 布 函数 。 


fti RRON uim Fn) uim Fc) , uz Fo) ,十 旱 特 征 值 的 
单 变 量 重 现 期 (7 表示 为 : 

. EY) . ED . EW 

! ]-F(x) ^? 1-F(x)  ? 1-F(x) 
式 中 并 为 干旱 事件 发 生 的 时 间 间 隔 ( 月 ); EQ) 为 干 
旱 事件 发 生 的 平均 时 间 间 隔 ( 月 )。 ED =N, HFP 
NN 为 干旱 的 时 间 序列 长 度 ( 月 );n 为 发 生 干 旱 事件 的 


(1) 


二 维 联合 重 现 期 可 表示 为 : 
T,= B 
i nx P(X, >x% UX, =x,) 
N (2) 


f = CLP. (s). (x) 
二 维 同 现 重 现 期 可 表示 为 : 


T = N 
and 

nx P(X, zx NX, zx) 
E N (3) 

nx{I-F,(x)-F,(x,)+ AOLA) 
三 维 联合 重 现 期 可 表示 为 : 
T= N 
n x P(X, 2x, UX zx, UX, =x) 
N (4) 


"ceo 
三 维 同 现 重 现 期 可 表示 为 : 


and 


nx P(X, >24 NX, zx, NX, =x) 


: 基于 CMIP6 模 式 的 叶 尔 羌 河 流域 未 来 水 文 干旱 风险 预 估 


N 
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nx fı Fu) - Fs) - F Qs) C[F s). Ea)  C[F. 2.6)  C[F. 9.9. (] - CLP. (x). (oe). F(x.) 


式 中 : 7, 为 干旱 要 素 的 联合 重 现 期 (a); 7 为 干 早 
要 素 的 同 现 重 现 期 (a) ; P(') 为 联合 分 布 函数 ; 
C[F. ()).F. (;)]  C[F, (x. F. (xa)] > CIE.. (2), F(x) 
为 二 维 联合 分 布 函 数 ; C[E Qs) EF (ns). (| 28 — 
维 联合 分 布 函 数 ; FG). F(a). FQ) 为 干旱 要 
素 的 边缘 分 布 函数 。 


2 结果 与 分 析 


HEC-HMS 模型 的 率 定 与 验证 
选取 2005 一 2010 年 的 叶 尔 羌 河流 域 出 山口 观 
测 数 据 对 HEC-HMS 模型 进行 率 定 ,并 利用 1986 一 
2004 年 和 2011 一 2014 年 2 个 时 期 的 数据 对 模型 进 
行 验证 (图 2)。 将 流域 模拟 的 日 径流 的 时 间 序 列 与 
实测 径流 的 时 间 序 列 进行 比较 ,以 NSE、R*、PBIAS 
等 参数 作为 评价 模拟 效果 的 指标 ( 表 3)。 从 图 2 和 
精度 验证 结果 中 可 以 看 出 ,该 模型 无 论 是 在 率 定 期 
还 是 在 验证 期 都 表现 良好 , 率 定 期 和 验证 期 的 NSE 
分 别 为 0.80.0.72 和 0.83 , 尼 值 分 别 为 0.81.0.70 和 
0.85 ,总 体 上 可 以 较 好 地 模拟 研究 区 的 径流 过 程 。 
2.2 叶 尔 芜 河流 域 未 来 流量 变化 趋势 预 估 

利用 HEC-HMS 模型 模拟 了 历史 与 未 来 3 个 
SSPs 情景 下 径流 变化 (图 3)。 结 果 表 明 ,未 来 3 种 


2.1 


3000 - 
~~~ SEMEL 证 
2500- 一 一 模拟 值 ii 
q, 2000 i 
$ 
ji 1500 | 
E 


1000 


500 


| 


(5) 

SSPs 情景 下 , 叶 尔 芜 河 径 流量 在 2015 一 2100 年 将 显 
著 增加 (P<0.1)。 在 不 同 SSPs 情景 下 ,未 来 径流 的 
模拟 结果 存在 一 定 的 变化 幅度 差异 ,但 所 有 情景 下 
的 径流 未 来 变化 趋势 一 致 。 具 体 而 言 ,在 SSP126、 
SSP245 和 SSP370 人 情景 下 叶 尔 芜 河流 域 流 量 将 分 别 
110.58 n?*s*a!,1.22 mesa ^ $12.21 m^: s "a! AY 
速率 增加 。 
2.3 时 尔 芜 河流 域 未 来 标准 化 径流 系数 的 时 间 变 
化 趋势 

基于 水 文 模型 HEC-HMS 模型 和 CMIP6 气 候 模 
式 数据 模拟 的 叶 尔 羌 河流 域 多 年 流量 ,计算 未 来 
(2015 一 2100 年 ) 不 同 SSPs 情 景 下 3 个 月 尺度 SRI 的 
时 间 变 化 (图 4)。 研 究 发 现 , 在 未 来 3 种 气候 情景 下 
2015$ 一 2100 年 叶 尔 芜 河流 域 上 游 的 SRI 均 呈 上 升 趋 
势 (P<0.1) ,上升 速 率 分 别 为 0.13.(10aj 0.21-(10a) ' , 
0.27…(10a)'。 男 外 ,未 来 3 种 气候 情景 下 , 叶 尔 羌 河 
流域 出 山口 SRI 的 最 小 值 将 分 别 出 现在 2037、2039 
年 和 2035 年 ,SRI 的 最 大 值 将 分 别 出 现在 2064、 
2080 年 和 2098 年 。 
2.4 流域 未 来 干旱 特征 变量 的 概率 分 布 

为 了 更 好 地 了 解 干旱 事件 的 总 体格 局 ,直观 地 
揭示 历史 和 未 来 水 文 干旱 事件 的 特征 差异 ,采用 小 
提琴 图 研究 未 来 (2015 一 2100 年 )4 个 干旱 特征 变量 
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图 2 IFK JET ERRE (2005—2010 4E ) IJ? EI (1 986—2004 ,2011—2014 4E ) SCHL AU 
Fig. 2 Observed and simulated averaged hydrographs for the calibration period (2005—2010) and validation 
period (1986—2004, 2011— 2014) of Yarkant River Basin 
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#3 HEC-HMS 模型 率 定 与 验证 精度 分 析 干旱 持续 时 间 和 干旱 强度 的 分 布 都 比较 分 散 ,存在 

Tab.3 Assessment of simulation results of the HEC-HMS 比较 明显 的 离散 值 ,干旱 烈度 和 烈度 峰值 的 分 布 比 
model during the calibration and validation periods 较 集 中 ;时 间 上 看 ,与 历史 时 期 (1961 一 2015 年 ) 相 
ME 年 份 NE R EE PAR SE TAD ARRE TEE RSPR AY 
率 定期 2005—2010 0.80 5.86 0.81 "T E di de mes 
验证 期 ”1986 一 2004 — 072 16.56 0.70 E e ur Es — P eae iy a 
2011—2014 0.83 9.74 0.85 下 未 来 时 期 的 干旱 烈度 较为 异常 ,这 主要 是 因为 这 

注 :NSE 为 纳什 系数 ;PBIAS 为 相对 偏差 ;为 决定 性 系数 。 一 时 期 存在 几 个 长 干旱 历时 且 干 时 峰值 较 大 的 干 


时 事件 。 干 旱 变量 的 变化 趋势 表明 ,从 历史 到 未 来 
于 时 持续 时 间 .干旱 强度 .干旱 烈度 和 干旱 烈度 峰 ”时 期 , 叶 尔 羌 河 流域 干旱 持续 时 间 .干旱 强度 和 列 
值 的 分 布 特征 。 从 图 5 中 可 以 看 出 , 叶 尔 芜 河流 域 。 上 度 峰 值 都 有 明显 上 升 ,干旱 程度 可 能 将 不 断 加 剧 。 
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注 :SSP 为 共享 社会 经 济 路 径 。 下 同 。 
图 3 1986 一 2100 年 叶 尔 羌 河 流域 年 径流 变化 趋势 
Fig.3 Variation trends of the annual runoff in Yarkant River Basin from 2015 to 2100 
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注 :SRI 为 标准 化 径流 系数 ;图 中 蓝 色 区 域 表 示 SRI>0 ,红色 区 域 表 示 SRI<0。 
图 4 未 来 SSP126.SSP245 .SSP370 情 景 下 叶 尔 羌 河流 域 SRI 变 化 趋势 
Fig.4 Variation trends of SRI in the Yarkant River Basin during future period under SSP126, SSP245 and SSP370 scenarios 
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图 5 历史 和 未 来 叶 尔 羌 河流 域 干旱 持续 时 间 干旱 强度 ,干旱 烈度 及 烈度 峰值 小 提琴 图 
Fig.5 Violin plots of the hydrological drought duration, drought severity, drought intensity and drought peak 


in Yarkant River Basin during historical and future periods 


2.5 未 来 多 变量 联合 干旱 概率 与 重 现 期 

单个 干旱 变量 的 分 析 不 足以 反映 流域 的 干旱 
特征 ,难以 对 未 来 水 文 干旱 发 生 的 概率 准确 预测 。 
为 了 精确 预 佑 未 来 叶 尔 羌 河 流域 上 游 水 文 干旱 发 
生 的 概率 , 评 佑 未 来 水 文 干旱 风险 ,本 研究 应 用 
Copula 函数 ,构建 水 文 于 旱 二 变量 .三 变量 联合 分 
布 ,对 未 来 不 同情 景 下 2015 一 2100 年 叶 尔 羌 河流 域 
水 文 干旱 的 联合 概率 进行 分 析 , 从 而 准确 地 反映 流 
域 未 来 的 水 文 干旱 和 情况。 根据 拟 合 效果 评估 结果 ， 
叶 尔 芜 河流 域 二 维 三维 Copula 联 合 概率 研究 选择 
干旱 持续 时 间 .干旱 烈度 和 烈度 峰值 3 个 变量 联合 
效果 最 好 ,因此 选择 这 3 个 变量 分 别 构建 相应 的 二 
HE 三维 Copula 联 合 分布 。 

根据 确定 的 最 佳 的 Copula 函数 形式 ,估算 了 相 
对 应 的 参数 ,构建 了 相应 的 二 维 三 维 联合 分 布 ( 表 
4)。 研 究 发 现 , 控 制 期 (1986 一 2014 年 ) 和 未 来 3 个 
SSPs 情 景 下 , 叶 尔 羌 河流 域 干旱 持续 时 间 和 干旱 烈 
度 的 联合 概率 分 别 为 0.37 .0.46 .0.35 .0.37; FEZ 
度 和 烈度 峰值 的 联合 概率 分 别 为 0.42 .0.43 ,0.43 , 
0.45 ;干旱 持续 时 间 和 烈度 峰值 的 联合 概率 分 别 为 


0.38 .0.41 .0.39 .0.40; 干 旱 持 续 时 间 干旱 烈度 和 和 烈 
度 峰 值 的 三 维 联合 概率 分 别 为 0.3 .0.37 .0.31 .0.35。 

为 探究 未 来 不 同情 景 下 叶 尔 羌 河 流域 水 文 十 
旱 概 率 变 化 趋势 ,将 未 来 不 同 SSPs 情 景 下 干旱 持续 
时 间 .干旱 烈度 和 烈度 峰值 的 三 维 联合 概率 进行 对 
比 (图 6)。 结 果 显 示 ,2015 一 2100 年 SSP126 ,SSP245 
和 SSP370 情景 下 水 文 干 旱 发 生 的 概率 都 呈 下 降 趋 
势 , 且 3 个 情景 下 水 文 干旱 事件 发 生 的 次 数 依 次 递 
减 ,干旱 概率 降低 的 速率 依次 递增 ,斜率 分 别 为 
-0.001 ,—0.007 ,-0.030.. 

对 识别 的 水 文 干旱 事件 进一步 进行 重 现 期 分 
析 ,以 研究 叶 尔 羌 河流 域 未 来 水 文 干旱 的 特征 。 给 
定单 变量 重 现 水 平分 别 为 2a、5 a、10 a、20a、50a 和 
100 a, 通 过 计算 单 变 量 干旱 特征 值 进而 得 到 两 变 
量 、 三 变量 联合 分 布 的 同 现 重 现 期 与 联合 重 现 期 
( 表 5)。 与 历史 时 期 重 现 期 相 比 , 叶 尔 羌 河流 域 未 
来 3 个 SSPs 情 景 下 ,不 同 单 变 量 重 现 期 所 对 应 的 烈 
度 峰 值 基本 都 低 于 历史 时 期 ,而 干旱 持续 时 间 基 本 
都 高 于 历史 时 期 。 以 单 变 量 重 现 期 为 10 a 为 例 , 历 
史 时 期 相对 应 的 干旱 持续 时 间 和 烈度 峰值 分 别 为 
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表 4 二 维 \ 三 维 Copula 函数 拟 合 参 数 
Tab.4 Fitting parameters of 2D and 3D Copula functions 
时 期 函数 参数 Dd&Di Di&Dp Dd&Dp Dd&Di&Dp 
控制 期 Gumbel 0 1.34 3.32 1.80 2.93* 
RMSE 0.12 0.09 0.11 0.05 
AIC —52.38 -61.21 -55.18 -342.51 
Frank 0 2.57* 13.42* 4.50* 9.41 
RMSE 0.12 0.08 0.11 0.05 
AIC -52.45 —64.16 —54.99 -334.57 
Clayton 0 0.30 3.50 0.53 3.85 
RMSE 0.13 0.09 0.12 0.06 
AIC —50.84 —61.34 -54.42 -323.45 
SSP126 Gumbel 0 2.59 2.15 1.57 3.51 
RMSE 0.05 0.04 0.06 0.09 
AIC —-347.79 -381.25 -333.49 -264.46 
Frank 0 11.14" 10.28* 4.80 12.08 
RMSE 0.04 0.03 0.04 0.10 
AIC —375.28 —408.83 -360.19 -261.79 
Clayton 0 0.05 0.08 1.94* 5.03* 
RMSE 0.14 0.11 0.04 0.09 
AIC —226.00 -253.63 -362.21 -271.12 
SSP245 Gumbel 0 1.14 2.72 1.29 2.23 
RMSE 0.04 0.04 0.05 0.11 
AIC —348.89 -334.98 -306.04 -231.66 
Frank 0 1.86" 11.99* 2.89 6.44" 
RMSE 0.03 0.03 0.04 0.11 
AIC -355.62 -349.88 -327.94 -232.54 
Clayton 0 0.02 0.07 1.32* 2.45 
RMSE 0.04 0.13 0.04 0.12 
AIC —-331.18 -213.75 -339.59 -221.57 
SSP370 Gumbel 0 1.46 4.15 1.75 3.00 
RMSE 0.07 0.05 0.07 0.06 
AIC -118.79 —138.45 —-123.35 —-130.76 
Frank 0 4.39* 14.43 6.44* 9.97* 
RMSE 0.06 0.05 0.05 0.06 
AIC —-126.57 —-136.16 —-134.79 —133.06 
Clayton 0 0.86 4.44* 1.16 4.00 
RMSE 0.07 0.04 0.07 0.06 
AIC —120.65 —-143.51 -119.68 -128.56 
注 :6 为 Copula 函 数 的 参数 ,RMSE 为 均 方 根 误差 ;AIC 为 赤 池 信息 准则 ; Dd&Di、Di&Dp、Dd&Dp 分 别 为 干旱 持续 时 间 和 干旱 烈度 \ 干 旱 烈度 和 


烈度 峰值 干旱 持续 时 间 和 烈度 峰值 的 二 维 联合 概率 ; Dd& Di& Dp 为 干旱 持续 时 间 、\ 干 旱 烈度 和 烈度 峰值 的 三 维 


合 效 果 最 好 。 


5.36 和 1.93; 未 来 时 期 在 SSP126 情 景 下 , 叶 尔 羌 河 
对 应 的 干旱 持续 时 间 和 烈度 峰值 分 别 为 589 和 
1.73;SSP245 情 景 下 ,对 应 的 干旱 持续 时 间 和 烈度 峰 
值 分 别 为 5.21 和 1.86;SSP370 情 景 下 ,对 应 的 干旱 


Ws ee N 
言 ,3 个 情景 下 干旱 持续 时 间 都 大 于 历史 时 期 ,而 列 

峰 q 换 句 话说 ,就 10 a 一 遇 的 
峰值 都 预计 将 低 于 历 


水 文 干旱 事件 而 言 , 未 来 烈度 
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图 6 叶 尔 羌 河 流域 未 来 3 个 SSPs 情 景 


三 维 联 合 干旱 概率 变化 趋势 


Fig.6 Variation trends of 3D joint probability of droughts in Yarkant River Basin during future period under three SSPs scenarios 


史 时 期 ,干旱 持续 时 间 预 计 将 高 于 历史 时 期 。 

根据 表 5 多 变量 联合 重 现 期 和 同 现 重 现 期 计算 
结果 ,结合 本 研究 水 文 干旱 等 级 划分 标准 可 知 : 叶 
尔 羌 河 未 来 3 个 SSPs 情景 下 烈度 峰值 最 小 分 别 为 
1.01 .0.9 .0.61 ,都 属于 轻 旱 范围 ,以 SSP126 Ay hil, 4 
单 变量 重 现 期 为 2a 时 , 叶 尔 羌 河流 域 Dd_Di、 
Di. Dp, Dd. Dp, Dd Di. Dp 的 同 现 重 现 期 分 别 为 
2.19 a,2.18 a、2.23 a,2.51 a, 基 本 一 致 ,可 见 发 生 持 
续 2.29 个 月 的 轻 旱 事件 的 重 现 期 约 为 2 a; 当 单 变量 
重 现 期 为 5a 时 , 叶 尔 羌 河 干旱 历时 约 为 4.12 个 月 ， 
烈度 峰值 约 为 1.59, 此 时 干旱 持续 时 间 和 烈度 峰值 
的 同 现 重 现 期 为 9.17 a, 可 见 叶 尔 羌 河 发 生 持 续 4.12 
个 月 的 中 旱 事件 的 重 现 期 约 为 9a, 与 历史 时 期 相 
比 ,发 生 中 旱 事 件 的 重 现 期 变 长 ;同样 地 , 单 变量 重 
现 期 为 10 a 时 ,历史 时 期 叶 尔 羌 河 发 生 持 续 5.36 个 
月 的 重 旱 事 件 的 重 现 期 约 为 19 a, 未 来 发 生 持 续 
5.21 个 月 的 重 旱 事件 的 重 现 期 约 为 28 a; 单 变量 重 
现 期 为 20 a 时 ,历史 时 期 叶 尔 羌 河 发 生 持 续 7.59 个 
月 的 极端 干旱 事件 的 重 现 期 约 60 a, 而 未 来 发 生 持 
续 6.97 个 月 的 极端 干旱 事件 的 重 现 期 约 115 a。 总 
体 而 言 , 与 历史 时 期 相 比 , 叶 尔 羌 河 未 来 各 类 干旱 
事件 的 重 现 期 都 将 有 不 同 程度 的 延长 。 


3 讨论 


影响 叶 尔 羌 河 流域 水 文 干 旱 的 因素 较 多 ,以 自 
然 因 素 为 主 ,如 温度 、 降 水 蒸发 和 冰雪 和 鹤 盖 ,人 类 


活动 对 上 游 水 文 干旱 的 影响 较 小 。 男 外 ,冰川 对 气 
候 变 化 的 响应 存在 滞后 效应 ,尤其 是 大 型 冰川 , 且 
对 于 冰川 年 际 变化 的 监测 也 存在 一 定 困难 。 考 虑 
到 冰川 的 积累 和 消融 也 主要 受降 雪 的 影响 , 故 通过 
对 融雪 的 监测 在 一 定 程度 上 可 以 来 表征 冰川 变化 。 
为 了 评价 自然 要 素 对 水 文 干旱 (径流 ) 的 相对 影响 ,本 
节选 择 了 4 个 主要 因子 , 即 从 Climatic Research Unit 
(CRU) 数 据 中 提取 1961 一 2015 年 的 年 温度 .降水 量 、 
Ag CRU 域 温 度 显著 升 高 ,降水 量 略 有 增加 。 温 度 升 
高 导致 山区 径流 增加 ,导致 融雪 和 冰川 融化 »。 叶 
尔 羌 河流 域 冰雪 融 水 对 河流 径流 的 贡献 在 1961 一 
2006 年 约 为 31.3% ,2000 年 后 增加 到 约 63.3% , 而 降 
水 补给 仅 占 22%25 。 在 以 冰雪 融化 为 主 的 流域 ,如 
叶 尔 羌 河流 域 , 径 流 受 温 度 的 影响 大 于 受降 水 的 
影响 。 

对 于 干旱 区 而 言 ,未 来 气候 变化 带 来 的 降水 增 
多 是 非常 有 限 的 ,对 叶 尔 羌 河流 域 未 来 3 个 SSPs 情 
景 的 气候 变化 预 估 结果 显示 ,年 降水 预计 将 从 历史 
时 期 的 83 mm 增加 到 未 来 的 92~107 mm, 相 对 增 量 
最 多 仅 20 mm ,温度 升 高 引起 的 冰川 积 雪 加 速 消 
融 对 水 资源 贡献 率 将 高 于 降水 对 水 资源 的 贡献 
率 。 对 预测 的 未 来 时 期 不 同 SSPs 情 景 下 温度 .降水 
和 径流 进行 相关 分 析 发 现 ,SSP126 情 景 下 径流 与 降 
水 相关 系数 (0.50) 高 于 与 温度 的 相关 系数 (0.49)， 
而 SSP245 和 SSP370 情 景 下 径流 与 温度 的 相关 系数 
(0.74) 高 于 与 降水 的 相关 系数 (0.63) ,这 一 关系 也 


806 tió WN 47 
表 5 叶 尔 芜 河流 域 历 史 与 未 来 时 期 重 现 期 
Tab.5 Return periods in Yarkant River Basin during historical and future periods 
时 期 干旱 变量 及 相应 重 现 其 
2a 5a 10a 20a 50a 100a 
历史 时 期 Dd/ 月 2.18 3.85 5.36 7.59 10.00 12.58 
Di 0.81 1.07 1.24 1.42 1.56 1.67 
Dp 0.91 1.43 1.93 2.58 3.19 3.77 
Dd DiT./a 2.09 5.91 14.02 37.17 82.25 161.24 
Dd e Ja 1.85 4.33 7.77 13.70 20.90 29.59 
Dp) Tana 241 8.33 23.89 76.49 190.98 405.30 
Di Dp)T./a 1.65 3.57 6.32 11.50 18.26 26.64 
Dd. Dp)T../a 2.27 7.32 19.82 60.32 146.20 304.69 
Dd Dp)T./a 1.72 .80 6.69 11.99 18.81 27.23 
Dd Di Dp)T../a 2.91 9.40 30.61 122.88 391.81 1000.00 
Dd Di Dp)T./a 1.68 3.56 5.92 9.84 14.64 20.41 
SSP126 Dd/H 2.29 4.12 5.21 6.97 10.23 18.48 
Di 0.78 1.24 1.45 1.62 1.74 1.83 
Dp 1.01 1.59 1.86 2.15 2.50 3.02 
Dd, Di)T../a 2.19 6.48 15.58 49.05 205.18 1562.50 
Dd Di)T./a 1.84 4.11 7.06 13.22 27.56 76.45 
Di. Dp)T../a 2.18 6.35 15.08 46.79 192.96 1452.99 
Di Dp)T./a 1.85 4.16 7.17 13.39 27.79 76.73 
Dd Dp)T../a 2.23 9.17 28.39 114.86 582.87 5054.03 
Dd Dp)T./a 1.81 3.46 5.86 11.45 25.35 73.95 
Dd Di Dp)T../a 2.51 9.20 24.57 82.12 346.75 2603.72 
Dd Di Dp)T./a 1.86 3.51 5.44 9.45 18.92 51.44 
SSP245 Dd/ 月 1.17 4.26 5.88 8.36 15.17 20.03 
Di 1.88 5.07 9.57 17.81 42.62 61.09 
Dp 0.90 1.40 1.73 2.08 2.63 2.88 
Dd Di)T../a 2.14 6.17 14.95 40.64 203.27 418.16 
Dd ih /a 1.88 4.31 7.67 13.51 31.19 44.92 
Dp) Panda 2.39 11.09 36.93 131.80 863.66 1909.51 
Di_Dp)T./a 1.73 3.29 5.88 10.99 27.92 41.45 
Dd. Dp)T../a 2.20 9.70 33.07 120.80 807.88 1794.67 
Dd Dp)T./a 1.84 3.43 5.99 11.07 27.98 41.50 
Dd Di Dp)T./a 2.40 17.66 169.82 594.86 1161.75 2129.31 
Dd, Di. Dp)T./a 1.74 3.41 5.63 9.46 21.15 30.27 
SSP370 Dd/ 月 1.85 3.77 9.84 22.98 69.41 132.10 
Di 1.10 3.53 10.73 25.50 74.90 138.81 
Dp 0.61 1.12 1.57 1.89 2.25 2.44 
Dd D1)T../a 3.67 527 12.97 36.15 175.58 484.10 
Dd ee Ja 3.65 4.71 8.08 13.74 29.96 49.36 
Dp) Panda 3.71 5.18 10.82 23.67 79.29 179.41 
Di Dp)T./a 3.61 4.79 921 17.18 37.80 59.70 
Dd. Dp)T../a 3.72 5.72 13.71 36.16 163.15 435.11 
Dd Dp)T./a 3.59 4.40 7.81 13.73 30.36 49.93 
Dd Di Dp)T,./a 3.71 5.42 13.95 40.72 206.81 567.91 
Dd Di Dp)T./a 3.59 4.35 7.65 12.87 25.64 39.51 
注 : Dd、Di 和 Dp 分 别 为 干旱 持续 时 间 、 干 旱 烈度 和 烈度 峰值 ;7 为 干旱 特征 变量 的 二 维 或 三 维 同 现 重 现 期 ;7, 为 二 维 或 三 维 联合 重 现 期 。 


在 其 他 干旱 区 流域 如 塔里木 河流 域 也 得 到 了 证 
明 "”。 叶 尔 羌 河 下 游 无 产 流 能 力 ,下 游 绿洲 生 态 系 
统 和 灌区 农业 用 水 全 部 依赖 于 上 游 山区 来 水 ,虽然 
生态 坝 的 建设 有 利于 生态 水 的 输送 ,可 以 在 一 定 程 
度 上 减少 干旱 灾害 的 威胁 "3, 但 水 文 干旱 比 其 他 自 
然 灾害 更 频繁 地 影响 着 下 游 及 整个 塔里木 河流 域 
的 生态 环境 和 人 类 活动 “”。 水 文 干旱 对 干旱 区 流 
域 生态 和 社会 经 济 发 展 的 威胁 最 大 ,对 水 文 干 旱 的 
研究 更 加 具有 现实 意义 。 

4 结论 

(1) HEC-HMS 模型 能 够 较 好 地 应 用 于 叶 尔 芜 
河流 域 ,在 未 来 3 个 SSPs 情景 下 ,2015 一 2100 年 叶 
尔 羌 河 流域 径流 将 显著 增加 (P<0.1), 叶 尔 羌 河 流域 
上 游 的 SRI 也 均 呈 显著 上 升 趋势 (P<0.1) ,上升 速 率 
分 别 为 0.13.(10aj .0.21.(10aj .0.27.(10aj-。 

(2) 就 水 文 干旱 特征 变量 而 言 , 与 历史 时 期 
(1961 一 2015 年 ) 相 比 , 叶 尔 羌 河 未 来 3 个 SSPs 情 景 
下 各 个 干旱 特征 变量 (持续 时 间 ,干旱 强度 .干旱 烈 
BE .烈度 峰值 ) 的 均值 和 中 位 数 基本 都 旦 上 升 趋势 ， 
未 来 干旱 事件 的 持续 时 间 .干旱 强度 .干旱 烈度 都 
将 增 大 。 

(3) 与 历史 时 期 (1961 一 2015 年 ) 相 比 ,未 来 3 个 
SSPs 情景 下 叶 尔 羌 河流 域 2015 一 2100 年 SSP126、 
SSP245 和 SSP370 情景 下 水 文 干 时 发 生 的 二 维 、 三 
维 联合 概率 都 将 呈 下 降 趋势 ,同时 干旱 事件 发 生 的 
重 现 期 部 将 有 不 同 程 度 的 延长 。 温 度 升 高 引起 的 
冰川 积 雪 加 速 消融 是 叶 尔 羌 河 上 游 来 水 量 增加 的 
主要 原因 ,虽然 总 体 看 来 未 来 径流 将 呈 增 多 趋势 ， 
水 文 干旱 概率 可 能 有 所 降低 ,但 气候 变化 带 来 的 不 
确定 性 增加 将 导致 水 文 干旱 事件 强度 等 特征 变量 
的 增 大 ,干旱 程度 可 能 较 历史 时 期 更 为 严重 。 
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Prediction of future hydrological drought risk in the Yarkant 
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Abstract: Global warming has led to the increased frequency of extreme events such as droughts, posing signifi- 
cant threats to ecological security and sustainable socioeconomic development, particularly in arid regions, which 
are highly sensitive and responsive to climate changes. This paper employs the distributed hydrological model 
HEC-HMS, utilizing observed meteorological and hydrological data from basin stations and global climate model 
data from the Sixth International Coupled Model Intercomparison Program (CMIP6), to simulate and forecast the 
historical (1986—2014) and future (2015—2100) runoff trends and hydrological drought risks in the Yarkant Riv- 
er Basin (an essential tributary of the Tarim River), Xinjiang, China. The findings indicate that: (1) The HEC- 
HMS model is well-suited for arid basin areas. Under the three shared socioeconomic pathways (SSPs) scenarios, 
the runoff and standardized runoff index (SRI) in the Yarkant River Basin are projected to significantly increase 
(P«0.1), with the SRI growth rate estimated at approximately 0.13-0.27 * (10a) '. (2) A comparative analysis of 
the marginal distributions of four drought characteristic variables in the basin for both historical and future peri- 
ods reveals that the duration and intensity of future droughts will exceed those in the historical record, with a con- 
tinuous rise in drought event magnitudes. (3) Moreover, the joint probability of future hydrological droughts in 
the Yarkant River Basin is expected to decrease relative to the historical period, leading to a prolonged return peri- 
od for future droughts. The outcomes of this study offer valuable scientific references for water resource manage- 
ment and the development of strategies to mitigate hydrological drought risks in the basin. 
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